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By means of a Monte Carlo cascade method, to account for the rescattering of the outgoing 
nucleon, we study the charged and neutral current inclusive one nucleon knockout reactions off 
nuclei induced by neutrinos. The nucleon emission process studied here is a clear signal for neutral- 
current neutrino driven reactions, and can be used in the analysis of future neutrino experiments. 
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I. INTRODUCTION AND THEORETICAL FRAMEWORK 

Recently an interest in neutrino scattering off nuclei has raised because of its implications in the experiments on 
neutrino oscillations based on large Cerenkov detectors. The presence of neutrinos, being chargeless particles, can 
only be inferred by detecting the secondary particles they create when colliding and interacting with matter. Nuclei 
are often used as neutrino detectors, thus a trustable interpretation of neutrino oscillation data heavily relies on 
detailed and quantitative knowledge of the features of the neutrino- nucleus interaction [lj. For instance, in the case of 
neutrino processes driven by the electroweak Neutral Current (NC), the energy spectrum and angular distribution of 
the ejected nucleons are the unique observables. Besides, the study of these distributions might also help to improve 
on our knowledge of the quark and gluon substructure of the nucleon, and in particular on the amount of nucleon 
spin carried by strange quarks . 

At intermediate energies, above the nuclear giant resonance and below the A(1232) regions, neutrino-nucleus 
interactions have been studied within several approaches. Several different Fermi gas, Random Phase Approximation 
(RPA) and shell model based calculations have been developed during the last 15 years 0-01 to compute neutrino 
or antineutrino induced single-nucleon emission cross sections. Most of the calculations use the plane wave and 
distorted wave impulse approximations (PWIA and DWIA, respectively), including or not relativistic effects. The 
PWIA constitutes a poor approximation, since it neglects all types of interactions between the ejected nucleon and the 
residual nuclear system. The DWIA describes the ejected nucleon as a solution of the Dirac or Schrodinger equation 
with an optical potential obtained by fitting elastic proton-nucleus scattering data. The imaginary part accounts for 
the absorption into unobserved channels. This scheme is incorrect to study nucleon emission processes where the state 
of the final nucleus is totally unobserved, and thus all final nuclear configurations, either in the discrete or on the 
continuum, contribute. The distortion of the nucleon wave function by a complex optical potential removes all events 
where the nucleons collide with other nucleons. Thus, in DWIA calculations, the nucleons that interact are lost when 
in the physical process they simply come off the nucleus with a different energy, angle, and maybe charge, and they 
should definitely be taken into account. A clear example which illustrates the deficiencies of the DWIA models is the 
neutron emission process: (vi,l~n). Within the impulse approximation neutrinos only interact via Charged Current 
(CC) interactions with neutrons and would emit protons, and therefore the DWIA will predict zero cross sections for 
CC one neutron knock-out reactions. However, the primary protons interact strongly with the medium and collide 
with other nucleons which are also ejected. As a consequence there is a reduction of the flux of high energy protons 
but a large number of secondary nucleons, many of them neutrons, of lower energies appear. 

In this talk, we present results for the QE (i/j,i/jJV), (yi,l~N), (vi,T>iN) and (Di,l + N) reactions in nuclei. We use 
a Monte Carlo (MC) simulation method to account for the rescattering of the outgoing nucleon. The first step is the 
gauge boson {W^ and Z° ) absorption in the nucleus, we take this reaction probability 1 from the microscopical many 
body framework developed in Refs. [l5llT^ | for CC and NC induced reactions. Starting from a Local Fermi Gas (LFG) 
picture of the nucleus, which automatically accounts for Pauli blocking, several nuclear effects are taken into account 
in that scheme: i) a correct energy balance, using the experimental Q— values, is enforced, ii) Coulomb distortion of 
the charged leptons is implemented by using the so called "modified effective momentum approximation" , Hi) medium 



1 It is given by the inclusive QE cross sections d 2 a/dQ'dE' (£7', E' are the solid angle and energy of the outgoing lepton, for a fixed 
incoming neutrino or antineutrino Laboratory (LAB) energy. We also compute differential cross sections with respect to d 3 r. Thus, we 
also know the point of the nucleus where the gauge boson was absorbed, and we can start from there our MC propagation of the ejected 
nucleon 
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polarization (RPA), including A— hole degrees of freedom and explicit pion and rho exchanges in the vector-isovector 
channel of the effective nucleon-nucleon force, and Short Range Correlation (SRC) effects arc computed, and finally 
iv) the nucleon propagators are dressed in the nuclear medium, which amounts to work with a LFG of interacting 
nucleons and it also accounts for reaction mechanisms where the gaug e boson, or Z° , is absorbed by two nucleons. 
This model is a natural extension of previous studies on electron |l7j|. photon |18| and pion |l9l l20j dynamics in nuclei 
and predicts QE neutrino-nucleus (differential and integrated) cross sections with an accuracy of about 10-15%, at 
intermediate energies |2"ill22| . 

After the absorption of the gauge boson, we follow the path of the ejected nucleon through its way out of the 
nucleus using a MC simulation method to account for the secondary collisions. Details on the MC simulation can be 
found in rial . This MC method was designed for single and multiple nucleon and pion emission reactions induced by 
pions |23t l24j and has been successfully em ploy ed to describe inclusive (7, 7r), (7, N), (7, NN),..., (7, Ntt),... [25ll26| . 
(e,e'7r), (e,e'N), (e, e'NN),..., (e, e'Nir),... |22| reactions in nuclei or the neutron and proton spectra from the decay 
of A hypernuclei (23. Thus, we are using a quite robust and well tested MC simulator. 



II. RESULTS AND CONCLUDING REMARKS 



The nucleon spectra produced by CC processes induced by muon neutrinos and antineutrinos of 500 MeV are shown 
in Fig.^for argon. Of course neutrinos only interact via CC with neutrons and would emit protons, but these primary 
protons interact strongly with the medium and collide with other nucleons which are also ejected. As a consequence 
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FIG. 1: Charged current 40 ' Ar(v^, pT +N) (upper panels) and 40 Ar(P M , /i + + N) (lower panels) cross sections as a function of the 
kinetic energy of the final nucleon for an incoming neutrino or antineutrino energy of 500 MeV. Left and right panels correspond 
to the emission of protons and neutrons respectively. The dashed histogram shows results without nucleon rescattering and 
the solid one the full model. 



there is a reduction of the flux of high energy protons but a large number of secondary nucleons, many of them 
neutrons, of lower energies appear. Our cascade model does not include the collisions of nucleons with kinetic energies 
below 30 MeV. Thus, the results at those low energies are not meaningful and are shown for illustrative purposes only 
in FigUl 
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The flux reduction due to the quasielastic NN interaction can be easily accommodated in optical potential calcu- 
lations. However in those calculations the nucleons that interact are lost when in the physical process they simply 
come off the nucleus with a different energy and angle, and may be charge, and they must be taken into account. 

The energy distributions of nucleons emitted after a NC interaction are shown in Figs. Inland In Fig. El we show 
the results for 40 Ar at two different energies. In both cases we find the large effect of the rescattering of the nucleons. 
For 500 MeV neutrinos the rescattering of the outgoing nucleon produces a depletion of the higher energies side of 
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FIG. 2: Neutral current 40 Ar(u, v + N) at 500 MeV (upper panels) and 150 MeV (lower panels) cross sections as a function of 
the kinetic energy of the final nucleon. Left and right panels correspond to the emission of protons and neutrons respectively. 
The dashed histogram shows results without rescattering and the solid one the full model. 



the spectrum, but the scattered nucleons clearly enhance the low energies region. For lower neutrino energies, most 
of the nucleons coming from nucleon nucleon collisions would show up at energies below the 30 MeV cut. 

As expected, the rescattering effect is smaller in lighter nuclei as can be seen in Fig. for oxygen. In all cases the 
final spectra of protons and neutrons are very similar. 

The ratio of proton to neutron QE cross section could be very sensitive to the strange quark axial form factor of 
the nucleon, and thus to the g 8 A parameter U S S [HI • The sensitivity to the collisions of the final nucleons is larger 
for both heavier nuclei and for larger energies of the neutrinos. In Fig. 0] one can clearly appreciate the importance 
of the secondary nucleons at the low energies side of the spectrum. 
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